WKB ANALYSIS FOR THE GROSS PITAEVSKII EQUATION 
WITH NON-TRIVIAL BOUNDARY CONDITIONS AT INFINITY 



THOMAS ALAZARD AND REMI CARLES 



Abstract. We consider the semi-classical limit for the Gross— Pitaevskii equa- 
tion. In order to consider non-trivial boundary conditions at infinity, we work 
in Zhidkov spaces rather than in Sobolev spaces. For the usual cubic nonlin- 
earity, we obtain a point-wise description of the wave function as the Planck 
constant goes to zero, so long as no singularity appears in the limit system. 
For a cubic-quintic nonlinearity, we show that working with analytic data may 
be necessary and sufficient to obtain a similar result. 



1. Introduction 

We study the semi-classical limit h — > for the Gross-Pitaevskii equation 

h 2 

ihdtu + —Au = Vu + f (\u\ 2 ) u, 
2m 

where x € W 1 . In the case of Bose-Einstein condensation (BEC), the external 
potential V = V(t, x) models an external trap, and the nonlinearity / describes the 
nonlinear interactions of the particles (see e.g. [10, 24, 18]). We consider two types 
of nonlinearity / (after renormalization) : 

• Cubic nonlinearity: f(\u\ 2 )u — (|u| 2 — l) u. 

• Cubic-quintic nonlinearity: f(\u\ 2 )u = (|u| 4 + A|m| 2 ) u, A e R. 

The cubic nonlinearity is certainly the most commonly used model in BEC. The 
defocusing nonlinearity corresponds to a positive scattering length, as in the case of 
87 Rb, 23 Na and 1 H. Note that this model is also used in superfluid theory. See e.g. 
[10, 24, 18] and references therein. The cubic-quintic nonlinearity, which is mostly 
used as an envelope equation in optics, is also considered in BEC for alkalimetal 
gases (see e.g. [13, 1, 23]), in which case A < 0. The cubic term corresponds to a 
negative scattering length, and the quintic term to a repulsive three-body elastic 
interaction. We also consider the case A > (positive scattering length). 



1.1. Cubic nonlinearity. Up to rescaling the Planck constant, we consider the 
limit e — > for: 

F 2 

(1.1) iedtu 6 + — Au e = Vu e + (\u e \ 2 -l)u e , x€R n ,n>l, 

(1.2) u £ (0,x) =a e { X y Mx)/e . 

Our initial data do not necessarily decay to zero at infinity. Typically, we do 
not assume a§ S L 2 (M. n ) (see Theorem 1.3 below). Recently, the Cauchy problem 
[11, 16] and the semi-classical limit [20] for (1.1) with V = have been studied more 
systematically. When the external potential V is zero, V = 0, the Hamiltonian 
structure yields, at least formally: 

j t (\\eVu%t)\\ 2 L2 + \\\u^t)\ 2 ~l\\ 2 L2 ) =0. 
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In this case, a natural space to study the Cauchy problem associated to (1.1) is the 
energy space (see e.g. [5, 16] and references therein) 

E = {uE ffio C (K") ; Vu G L 2 (R"), \u\ 2 - 1 G L 2 (R")}. 

For this quantity to be well defined, one cannot assume that u £ is in L 2 (R™); 
morally, the modulus of u £ goes to one at infinity. To study solutions which are 
bounded, but not in L 2 (R"), P. E. Zhidkov introduced in the one-dimensional case 
in [28] (see also [29]): 

(1.3) X"(R n ) = {u G Z,°°(R") ; Vu G i? 8 " 1 ^")}, s > n/2. 

We also denote X°° := n s> „/ 2 A s . The study of these spaces was generalized in the 
multidimensional case by C. Gallo [11]. They make it possible to consider solutions 
to (1.1) whose modulus has a non-zero limit as \x\ — > oo, but not necessarily 
satisfying \u £ (t, -)| 2 — 1 G L 2 (M. n ). We shall also use these spaces. 

Recently, P. Gerard [16] has solved the Cauchy problem for the Gross-Pitaevskii 
equation in the more natural space E, in space dimensions two and three. The 
main novelty consists in working with distances instead of norms, in order to apply 
a fixed point argument in E. In particular, the constraint \u £ (t, -)| 2 — 1 G L 2 (W n ) 
is satisfied. 

To our knowledge, if the initial data do not vanish at infinity, the introduction 
of an (unbounded) external potential in Gross-Pitaevskii equation has no phys- 
ical motivation. Note also that if V is an harmonic potential, then the formal 
Hamiltonian corresponding to (1.1) is necessarily infinite (see § 2.3). On the other 
hand, introducing a quadratic external potential or considering a quadratic initial 
phase (j>o makes no difference in our analysis. The model (1.1)— (1.2) with V = 
and 4>o quadratic is certainly more physically relevant, and does not seem to enter 
into the framework of the previous mathematical studies. Another motivation to 
introduce this external potential stems from the study of the semi-classical limit 
of the Schrodinger-Poisson system, where |u e | 2 — 1 is replaced with V £ given by 
AV £ = q (\u £ \ 2 — c). This models appears in the semi-conductor theory where the 
real number q models a charge, which we may take equal to one here, and the func- 
tion c = c(x) models a doping profile, which we may take to be c = 1. As in [2], we 
will prove that if V grows quadratically in space, then if \u £ (t — 0, -)| 2 — 1 G L 2 (R n ), 
one must not expect \u £ (t, - )| 2 - 1 G L 2 (R") for t > 0. 

Assumptions. We assume that the potential and the initial phase are of the form: 

• V e C°°(R t x R"), and V = V quad + V iin , where V quad {t,x) = t xM{t)x is 
a quadratic form, with M{t) G <S„(R) a symmetric nx n matrix, depending 
smoothly on t, and VVjin € C°°(Ru X s ) for all s > n/2. 

• 0o € C°°(R") 7 and </> = </> qua d + <?W, where 4> qua d(x) = t xQ n x is a qua- 
dratic form, with Qo a symmetric matrix in M nX n(M), and V0i; n G X°° . 

Note that our assumptions include the case where Vu n and </>n n are linear in x. 
In general, these functions are sub-linear in x, since their gradient is bounded. 

Lemma 1.1. There exist T > and a unique solution (j)^ G C°°([0, T] x R") to: 

(1.4) <9 t cik + ^IV^oikl 2 + V^uad = | c ik|t=O = 0quad • 

Moreover, ^> e ik is a quadratic form in x: 

(1.5) (j>eik(t,x) = *xQ(t)x, 
where Q(t) G <S n (R) is a smooth function oft. 
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Proof. Existence and uniqueness follow from [8, Lemma 1]. To prove that </> c ik is 
quadratic in x, seek e ik of the form (1.5). Then (1.4) is equivalent to the system 
of ordinary differential equations 

Q(t) + 2Q(t) 2 + M(t) = ; Q(0) = Q . 

The lemma then follows from Cauchy-Lipschitz Theorem. □ 

Remark 1.2. As in [2], we shall use the following geometrical interpretation of the 
above lemma. The time T is such that for t e [0, T], the map given by 

d t x(t,y) = V x <p eik (t,x(t,y)) = Q(t)x(t,y) ; x(0,y)=y, 

defines a global diffcomorphism on R™. Therefore, the characteristics associated to 
the operator d t + V(/> c ik • V do not meet for t € [0, T], and this operator is a smooth 
transport operator: 

(d t + V0 cik • V) (f(t, y)) = d t f (t, x(t, y)) . 
Note that if Q(t) and its anti-derivative commute, then we have 



x(t,y) = cxp 



/ Q(r)dn 
Jo 



Theorem 1.3. Suppose that there exist ao,ai e X°° such that: 
(1.6) ||a§ — a — £ai|| xs = o(e), Vs > n/2. 

There exist T* G]0,T] independent of e <=]0, 1], and a unique solution u e G C°° Ci 
L°°([0,T*] x R") to (1.1)-(1.2). Moreover, there exist a, (f> £ C°°([0,T*] x R") with 
a,V(j) € C([0,T*];X s ) for all s > n/2, such that: 



(1.7) limsup 



u e (t, ■) — a(i, -)e 



i(<j>(t,-)+<j> eik (t,-))/e 



L°°(R") 



0(t) ast^O. 



The functions a and (f> depend nonlinearly on 4>o and oq (see (3.1) below). There 
exists <j>W G L°°([0,T*] x R"), real-valued, with V4> (1) G C([0, T*]; X s ) for all 
s > n/2, such that: 



(1.8) limsup u^-ae^ e i( * + * cik)/£ =0. 

e^O L~([0,T,]xl») 

The modulation (p^ is a nonlinear function of </>o, &o and a\ (see (3.2) below). 
Remark 1.4. Several applications of this general results are given, in §3, §4 and §5. 
Remark 1.5. If we assume moreover 

\\a £ - a - sai\\ Xs = 0(e 2 ), Vs > n/2, 
then the above error estimate can be improved: 

0(e). 



ae i0 (1) e i(0+0cik)/£ 



L°°([0,T«]xR") 



Remark 1.6. The above result and System (3.2) below show that in general, it is 
necessary to know the initial amplitude a§ up to the order o(e) to describe the 
leading order behavior of the wave function u £ . It is not necessary to know 
with such precision to study the convergence of quadratic observables. See §6. In 
particular, in Theorem 6.1, we extend the result of [20] to the three-dimensional 
case (on a bounded domain, or outside a bounded domain). 
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Remark 1.7. Most of the results that we present here remain valid in a space- 
periodic setting, that is if we assume i£T". In that case, compactness arguments 
show that the proof of Theorem 1.3 remains valid when V G C°°(Rt x T") and 
<f>o € C°°(T"). On the other hand, the discussions in §2.3 and §5 become irrelevant 
on the torus. Finally, note that it is equivalent to work in Sobolev spaces, since 
X s {T n ) = H s (T n ) for s > n/2. 

The analysis detailed in §2 and §3 shows that the formal part of [7] can be 
justified in the present framework. We shall only state a typical consequence of 
this approach: 

Corollary 1.8 (Instability). Let n ^ 1, o ,ai G C^nX 00 ^"), withRe(a ai) 0, 
and <f)o € C°°(IR n ), with V<fio G X°° . Let u £ and v £ solve the initial value problems: 



ied t u e + -Au £ = (\u £ \ 2 - 1) u £ ; u £ \ t=0 = a o e J0o/e . 
e 2 

ied t v £ + -Av e = (\v £ \ 2 - 1) v e ; v £ \ t=Q = (a Q + S £ ai ) e* 



where S £ — > 0. ^ssitme i/iat i/iere ermfc ]VeN such that S £ /e 1 « — ► +oo. T/ien 
we caw find t £ — > smc/i i/iaf liminf ||u e (i e ) — ^ e (i e )||£cc > 0. In particular, 



.. . r ll U W llL°°([0.t e ]xR") 

lirnini 11 ; 1 = +oo. 



«|t=0 v \t=o LOO(R „ ) 
Remark 1.9. Note that if </> = 0, then we also have: 



lim inf ■ 



lL°°([0,t £ ]xR«) 



u \t=0 V \t=0 



+oo, Vs > n/2. 



This shows that the instability mechanism is not due to regularity issues. It is due 
to the fact that (1.1) is super-critical as far as WKB analysis is concerned: the 
small initial perturbation (of order 5 s ) yields a high-frequency perturbation of the 
evolution (a multiplicative factor of the form e~ 2lts Rc Oo a i)/ £ ). 

1.2. Cubic-quintic nonlinearity. Denote f\(y) = y 2 + Xy. We now consider 

(1 | isd t u £ + ^Au £ = f\ (\u £ \ 2 ) u £ , x£l",n)l, 

t u £ (0,x) = a £ o (x)e i4 ' o( - x ^ £ . 

Note that in (1.9), we assume that there is no external potential, V = 0. We 
also assume that there is no initial quadratic oscillation: <fi G C°°(IR ra ; R), with 
V0 O G The case A > 0, V ^ 0, with a§ G H°°, is contained in [8]. We assume 

^quad = here in order to consider non-zero boundary conditions at infinity. We 
also assume Vii„ = for simplicity only. 

Plugging an approximate solution of the form u £ « ae l< ^^ £ , with a and <f> inde- 
pendent of e, and passing to the limit e — > as in [14, 20], we find formally that 
(p,v) := (|a| 2 ,V</>) solves: 



(1.10) 



d t p + div (pv) = 0. 

d t v + v Vv + V(f x (p)) =0. 



If A > 0, then the problem is hyperbolic. Essentially, the result of Theorem 1.3 
remains valid. When A < 0, the above problem is hyperbolic for p > |A|/2 and 
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elliptic for p < |A|/2. This feature is reminiscent of Eulcr equations of gas dynamics 
in Lagrangian coordinates: 



(1.11) 



d t u + d x v = 0. 
d t v + 8 X (p(u)) = 0. 



As recalled in [22] , a typical mathematical example for van dcr Waals state laws is 
given by p(u) = (u 2 — l)u. The problem is hyperbolic if u > 1/V3, and elliptic if 
u < Hadamard's argument implies that the only reasonable framework to 

study (1.10) or (1.11) is that of analytic functions (see [22]). In this case, we refer 
to the approach of [15, 27]. More details are given in §7. When the elliptic region 
for (1.10) is avoided, then essentially, Theorem 1.3 remains valid: 

Theorem 1.10. Suppose that there exist ao,a\ € X°° such that: 

\\a £ Q - a - eai\\ xe = o(e), Vs > n/2. 



Assume moreover that 
• Either A > 0, 



with V0o G X°°, and: 



Or A < and there exists S > such that |a (a;)| 2 ^ S + Vx G 



Then there exist e*,T* > 0, and a unique solution u e 
to (1.9) /or a/Z e g]0, £*]. Moreover, there exist a,<j 
a,V<f> e C([0,T*];X s ) /or a/Z s > n/2, sucA t/iat; 



€ C 00 n L°°([0,T*] x M") 
G C°°([0,T*] x R") with 



lim sup 



u e (t, •) — a(t, -)e 



■■4{t,-)/e 



= 0{t) ast^O. 



There exists 4>^ G L°°([0,T,] x W 1 ), real-valued, with V<t>^ g C([0, T„]; X s ) /or 
all s > n/2, such that: 



lim sup 

e^0 



ae 



L~([0,T»]xM™) 



0. 



1.3. Structure of the paper. In §2, we construct the solution u e as u e 



a e 



j 

where a £ is complex-valued and $ e is real- valued. This yields the existence part of 
Theorems 1.3 and 1.10. The proof of these theorems is completed in §3, where the 
limit of (a e , <I> e ) as £ goes to zero is studied. We give three examples of applications 
of Theorem 1.3 in §4, in the case ^> e ik = 0. In §5, we study the time evolution of 
a non-trivial boundary condition at infinity when ^> e ik 7^ 0. In §6, we investigate 
the limit of the position and current densities. Finally, we explain why working 
in an analytic setting is often necessary (and always sufficient) in the case of the 
cubic-quintic nonlinearity. 



2. Construction of the solution 

2.1. Phase-amplitude representation: the case (f)^ = V = 0. When V and 
4>o are identically zero, the existence and uniqueness part of Theorem 1.3 was 
established by C. Gallo [11]. Note however that with our scaling, the fact that T* 
is independent of e g]0, 1] does not follow from [11]. Since the approach in Zhidkov 
spaces is rather similar to the one in Sobolev spaces, we shall essentially explain 
the new aspects of the proof. To treat both cubic and cubic-quintic nonlinearities, 
consider the general equation 

f(\u E \ 2 )u E , x g R™, n > 1, 
a £ (xy Mx)/£ , 



(2.1) 



ied t u £ + y Au e = 
u £ {0,x) = 
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where / £ C° 



). We keep the hierarchy introduced by E. Grenier [17]: seek 



u e = a € e , where a 6 is complex- valued, and $ e is real-valued. We impose 



(2.2) 



I|V$ £ | 2 



o 



t=o 



d, 



-a e A$ e 
2 



«-Aa 
2 



it=o 



As an intermediary unknown function, introduce the "velocity" v £ — V<I> e . Separate 
real and imaginary parts of a e , a £ = a\+ iaf,, and introduce: 

/ of \ / Re(a§) \ 



u = 



Un 



V < J 



Im(ag) 

V d n4>0 ) 



-A ... 
L = [ A ... 

0„ xn 



and 4(11,0 = $ZAKu)6 = 

i=i 



v-Z 
v£ f-^ 



«i tt 

2 ^ 

a 2 t 



where /' stands for /'(|ai| 2 + |a.2| 2 )- We now have the system: 
(2.3) 



a t u E + ^A J (u £ )^ = -Lu £ ; uf t=0 = u§. 



The matrices A,- are symmetrized by the matrix 



h 



which is symmetric positive if and only if /' (|ai| 2 + |a 2 | 2 ) > 0: this includes the 
case of the decofusing cubic nonlincarity (1.1), of the cubic-quintic nonlinearity 
(1.9) with A > 0, and of the cubic-quintic nonlincarity (1.9) with A < 0, provided 
that |ai| 2 + |a 2 | 2 > |A|/2. 

Proposition 2.1. Assume that Uq is bounded in X s for all s > n/2, uniformly for 
e £ [0, 1], and that there exists e* > and 5 > such that 

/'(Kl 2 ) > S > 0, Vi6R n , Vee[0,4 

Then for s > n/2 + 2, there exist T* > and a unique solution u £ £ C([0, T*]; X s ) 
to (2.3) for all e £ [0, £*]. In addition, this solution is in C([0, T*]; X m ) for all 
m > n/2, with bounds independent of e £ [0, £*]. 

Proof. Let s > n/2 + 2. As usual, the main point consists in obtaining a priori 
estimates for the system (2.3), so we shall focus our attention on this aspect. We 
have an a priori bound for u £ in L°° : 

l|u c (t)||L~« i ~+ / Y,\\ A ^)d J n £ (r)\\ L ^dr+ / || Au £ (t)|| l ^t 
Jo J=l Jo 

<K|| i0 o + f F{\\u e {T)\\ L oo)\\Vu e {T)\\ H .- 1 dT 

Jo 

ft 



+ C 



We infer: 
(2.4) 



/ ||Au e (r)|| ff3 - 2 dr. 
Jo 

U%t)\\ L ~ ^\K\\ L ~+ fGiWu^T^Wu^T^xsdT. 

Jo 
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To have a closed system of estimates, introduce P = (I — A)( S_1 )/ 2 V, so that 
||/lk.»||/||L~ + ||P/|| L2 . Denote 



if, 9) = / f(x)g(x)dx, 

the scalar product in L 2 . Since S is symmetric, we have 

^-(SPu £ (t),Pu £ (t)) = (d t SPu £ (t),Pu £ (t))+2Rc(Sd t Pu £ (t),Pu £ (t)), 
at 

So long as 

(2-5) /'(l« £ | 2 )^>0, 
we have the following set of estimates. First, 

(d t SPu £ (t),Pu £ (t)} < \\d t S\\ L ~\\Pu%t)\\ 2 L2 

<C 4 (||u e (t)|Ucc)||ftu e (t)|| L oc||u e (t)||^. 

Directly from (2.3), we have: 

\\d t u £ (t)\\ L ~ < C(||u £ (i)|k=o) ||Vu e (t)|Uoo + ||Au £ (t)|| L =o 

<C(\\u £ (t)\\ X s)\\u £ (t)\\ X s. 

Since SL is skew-symmetric, we have 

Rc{SLPu £ {t),Pu £ {t)) = 0, 

which prevents any loss of regularity in the estimates. For the quasi-linear term 
involving the matrices Aj, we note that since SAj is symmetric, commutator esti- 
mates (see [19]) yield: 

n 

£(SP(A,(u £ )dX),Pu £ (t)) ^C(\\u £ (t)\\ L ~)\\Pu £ (t)\\U\Vu £ (t)\\ L ~ 

^C(\\u £ (t)\\ Xs )\\Pu £ (t)\\ 2 L2 . 

Finally, we have: 

|(W(t),Pu £ (t)) ^C(\\u £ (t)\\ Xs )\\u £ (t)\\ 2 xs . 

This estimate, along with (2.4), shows that on a sufficiently small time interval 
[0, T»], with T„ > independent of e G [0, (2.5) holds. This yields the first part 
of Proposition 2.1. 

The fact that the local existence time docs not depend on s > n/2 + 2 follows 
from the continuation principle based on Moser's calculus and tame estimates (see 
e.g. [21, Section 2.2] or [26, Section 16.1]). □ 

The existence part of Theorem 1.10 and of Theorem 1.3 when <j> e ^ = follows. 
Indeed, define $ £ by 

= <h - f* QK(T)| 2 + / (K(t)| 2 )) dr. 

We check that <9 t (V$ £ - v £ ) = V<9 t $ £ - d t v £ = 0, so that V$ £ = v £ , and ($ e , a £ ) 
solves (2.2). Finally, uniqueness for (2.1) follows from energy estimates. If u £ ,v £ £ 
C°° n L°°([0, T»] x R") solve (2.1), then w £ := u £ - v £ satisfies: 

ied t w £ + £ ^Aw £ = f (\u £ \ 2 ) u £ ~f (\v £ \ 2 ) v £ ; wf t=0 = 0. 
We have, for t G [0,T„], 

||w £ ||l,oo(o it; L2) C (||u £ ||l~([0,T,]xM"), II w£ |U°°([0,T»]xR")) ||w £ ||Li(0,t;L 2 ), 
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and Gronwall lemma yields w £ = 0. 



2.2. Phase-amplitude representation: the case o ik ^ 0. We know consider 
(1.1)— (1.2) only: the nonlinearity is exactly cubic. To take the presence of V 
and </i>quad into account, we proceed as in [8]: we construct the solution as u £ — 
a e e i(0 e +^eik)/e. analogue of (2.2) is: 



d t $ e + - |V$ £ | 2 + V + \a £ \ 2 -1 = 



<j) e 



t=o 



d t a £ + V$ £ • Va £ + -a e A$ e 



i-Aa £ 
2 



\t=o 



Set 4> e = (f> £ + c ik- The introduction of ^ e ik allows us to get rid of the terms 
Vq Ua d and f^quad, and work in Zhidkov spaces. The above problem reads, in terms 
of(0 e ,a £ ): 



(2.6) 



+ \ IWf + V0 oik • V4> £ + V lin + \a £ \ 2 



d t a £ + \7cj) £ ■ Va £ + V</> ci k ■ Vc 



ia e A0 £ + ia £ A0 cik 



0. 



i^Aa £ , 



<t> £ 



= 9hn 



lt=0 



a. 



Resume the notations of the previous paragraph, with now: 
I ° \ 

/ v£ 

S= d ^ , and 4(u,0 ( 

: i=i 

The system (2.3) is replaced by: 







o ■ 



<ji tc 
2 S 
a 2 t 



2ai £ 2a 2 £ w • £/„ 



(2.7) <9 t u £ + ^ A^u^u 5 + V^ik • Vu E + Mu £ + S = -£u e 



u |t=o 



u 



0' 



where M = M(t) is a matrix depending on time only, since 4> e ^ is exactly quadratic 
in x. This aspect seems necessary in the proof of Proposition 2.2 below. This 
explains why we make Assumptions 1.1, and do not content ourselves with general 
sub-quadratic potential and initial phase as in [8]. The important aspect to notice 
is that since the nonlinearity in (1.1) is exactly cubic, then the matrices Aj are 
symmetrized by a constant matrix, namely: 



S = 



h 






It 

1 ^ // 



In [8], nonlinearities which are cubic at the origin were considered (as in [17]), and 
the possibly quadratic phase e ik made the assumption xoq <E L 2 (R n ) apparently 
necessary, to control the time derivative of the symmetrizer. Of course, we want to 
avoid this decay assumption for the Gross-Pitaevskii equation, so working with a 
constant symmetrizer is important. 

Proposition 2.2. Assume that ug is bounded in X s for all s > n/2, uniformly for 
s € [0, 1]. Then for s > n/2 + 2, there exist T* €]0, T], independent of e G [0, 1], 
and a unique solution u £ € C([0, T*]; X s ) to (2.7). In addition, this solution is in 
C([0, T*};X m ) for all m > n/2, with bounds independent of e € [0, 1]. 

Sketch of the proof. The proof follows the same lines as the proof of Proposition 2.1, 
so we shall only point out the differences. 
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Let s > n/2 + 2. By construction, the operator d t + V</> ; k • V is a transport 
operator along the characteristics associated to (/> c i k , which do not intersect for 
t e [0, T]. Therefore, we have an a priori bound for u e in L°°: 

||u e (t)|| L ^ ^|| U g|| L oo + / ^T\\A 3 {u)d 3 u{T)\\ L ~dT 

Jo j=1 

+ J (c||u £ (r)|| L oo + ||S(r)|| L oo + ||Au £ (r)|| L oc)rfT 

(2.8) + C f (1 + ||u e (r)||xO \\^(r)\\ xs dr + ||S|| L ^ ([0iT];Xs) . 

Jo 

To have a closed system of estimates, resume the operator P = (I - A)( 5 " 1 )/ 2 V, 
so that \\f\\ X s » ||/|| L oo + ||P/|| L 2. We have 

j t (SPu £ (t), Pu £ (t)) = 2Rc (Sd t Pu £ (t), Pu £ (t)) , 

since S is constant symmetric. Since SL is skew-symmetric, we have 

Rc(SLPu £ (t),Pu £ (t)) = 0. 

For the quasi-linear term involving the matrices Aj, we note that since SAj is 
symmetric, commutator estimates yield: 

n 

J=l 

Next, write 

(SP (V0 oik • Vu £ (i)) , Fu £ (t)) = (5V0 cik • VPu £ (i), Pu £ (t)) 

+ (S[P,V0 oik -V]u e (t),Pu £ (t)>. 

The first term of the right-hand side is estimated thanks to an integration by parts: 
2 Re (5V0 cik • VPu £ (t), Pu £ (t)) = J SV</> cik (t, x) • V|Pu £ (t, x)\ 2 dx 

SA(j) cik (t, x)\Pu £ (t, x)\ 2 dx. 



For the second term, we notice that [P, V0 o i k • V] = ipV, where ip = ip(t,D) is a 
pseudo-differential operator in x, of order s — 1, depending smoothly of t e [0, T]. 
Therefore, 

2 Re (SP (V0 cik • Vu £ (t)) , Pu £ (t)) < ||u e (i)|&. . 

The fact that M is independent of x is crucial here, to ensure that P(Mu e ) e L 2 
for u £ G X s . If ikf depended on x, that is if (/> c ; k was not a polynomial of order at 
most two, the low frequencies might be a problem at this step of the proof. Finally, 
we have: 

|<W(t),Pu e (t)) <C(||u £ (t)|| Xs )IIK(t)||^. 

This estimate, along with (2.8), yields the first part of Proposition 2.2. We conclude 
like in the proof of Proposition 2.1. □ 

The existence part of Theorem 1.3 follows from the above result, by setting 
<p £ (t) = fa ~ I (\ KMI 2 + Wcik(r) • v £ (t) + Hi»(T) + |a £ (r)| 2 - dr. 
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Finally, uniqueness for (1.1)— (1.2) follows from energy estimates. If u £ ,v £ G C°° n 
L°°([0,T*] x R n ) solve (1.1) (1.2), then w E := u E - v £ satisfies: 

e 2 

ied t w £ + yAw £ = (V - l)w £ + \u e \ 2 u e - \v £ \ 2 v £ ; tuf t=0 = 0. 
We have, for t G [0,T«], 

l|w £ ||L=o (0it . L 2 ) < (iKHl^Qo.T.lxR") + II w£ Hl-([0,T»]xR")) II^IU^O,*;!, 2 ), 

and Gronwall lemma yields w e = 0. 

2.3. On the Hamiltonian structure. When V — V(x) is time-independent, 
(1.1) formally has a Hamiltonian structure, with 

H=h\eVu e (t)\\h+ [ V(x)\u%t,x)\ 2 dx+\\\\u%t)\ 2 -l\\l 2 . 

When V = 0, this structure is used in [16] to prove the global existence of solutions 
in the energy space. On the other hand, suppose that V is, say, harmonic: 

n 

J=l 

where the constants Xj are not all equal to zero. Then necessarily, H is infinite: 
suppose for instance that \i > 0. Then if d xi u e (t, -),X\U e (t, •) G L 2 (R n ), the un- 
certainty principle (a simple integration by parts, plus Cauchy-Schwarz inequality 
in this case) yields: 

u e {t,-) G L 2 (R n ). 

Therefore, the constraint \u £ (t, -)| 2 — 1 G L 2 (R n ) cannot be satisfied, for otherwise, 
1 = 1- K(t, -)| 2 + \u e {t, -)| 2 G L 2 (R n ) + L^W 1 ). 

Similarly, assume that V = 0, but </> qua d 7^ 0: rapid quadratic oscillations arc 
present in the initial data. We have 

eWu\ t=0 = (eWa e + ia e V<p ) e^°/ £ . 

Therefore, the above quantity is in L 2 provided that Vag, agV^quad G L 2 (R n ). If for 
instance (j)q Ua <i(x) — cx\ with c ^ 0, the last assumption means that x\a\ G L 2 (R n ), 
which brings us back to the previous discussion. 

We shall see in Section 5 that if ^ e ik ^ 0, and if ag G X°° is such that 

|a e | 2 - 1 G L 2 (K"), 

then the last constraint present in H is not propagated in general. In small time 
at least, one has generically 

\u e {t,-)\ 2 -l^L 2 (R n ). 

3. Semi-classical analysis 

We now complete the proof of Theorem 1.3. The end of the proof of Theorem 1.10 
follows essentially the same lines, so we omit it. The main adaptation is due to 
the fact that when the nonlincarity is not exactly cubic, the symmetrizer S is not 
constant. We refer to [17] or [8], to see that the proof below is easily adapted. 

Introduce (<j),a), solution to (2.6) with e = 0, that is 

dt<f> + \ |V0| 2 + V0 oik • V0 + Vu n + \a\ 2 -1 = 0; <f>\ = 0i in , 
(3-1) \ ! i 

d t a + • Va + V0 oik • Va + -aA<f> + -aA^ cik = ; o| t=Q = a . 
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It is a particular case of Proposition 2.2 that (3.1) has a unique solution, such that 
a,V0G C([0,T*];X s ) for all s > n/2. 

Proposition 3.1. Under the assumptions of Theorem 1.3, let (<j) 6 ,a £ ) and ((f>,a) 
be given by (2.6) and (3.1) respectively. For all s > n/2, there exists C s such that 

l|V(</> £ - 0)||l-([o,t»]:A=) + \\a £ - a|U~([o,T,] ; x-) ^ C s e. 

Sketch of the proof. We shall give the outline of the proof, since it is very similar 
to the case of Sobolev spaces [8] . The differences are those pointed out in the proof 
of Proposition 2.2. Resuming the notations of §2, set 



u 



/ Re a \ 
Ima 

dx<f> 



/ Re(a^ - a ) \ 
Im(a§ - a ) 




V d n <P J 
Denoting w £ = u e u, (2.7) yields: 

n 

d t w E + V (A,(u e )a,u e - A 3 (u)d 3 u) + V^ cik • Vw E + Mw e = ^Lw £ + ^-Lu , 

3=1 

wf t=0 = w§. 

We know by Proposition 2.2 that u e and u are bounded in C([0, T*]; X s ) for all 
s > n/2. The source term £ in (2.7) is now replaced by §£u, which is of order 
O(e) in C([0,T*];X s ), and we have easily, for s > n/2 and t G [0,T*]: 

l|w £ |U=o ([0;t];XS) \\w s \\ X s + 0(e) + [ \\^(r)\\ xs dr. 

Jo 

The proposition follows from Gronwall lemma. □ 

Remark 3.2. Note that for the time T* in Proposition 3.1 (as well as in Proposi- 
tion 3.4 below), we can pick the life-span of (<p, a), the solution of (3.1). Indeed, 
the error estimate and the standard continuity argument show that (<fi £ , a £ ) cannot 
blow-up as long as (<p, a) remains smooth, provided that e is chosen sufficiently 
small. In particular, if ((f), a) remains smooth globally in time, then for any r > 0, 
we can find e(r) > such that Proposition 3.1 and Proposition 3.4 below remain 
valid on [0,t] for e <E]0,e(r)]. On the other hand, one must not expect T» = oo in 
general: the solution to (6.1) may not remain smooth for all time. See [25]. 

Corollary 3.3. There exists C such that for all t e [0, T*], 

||0 e (t,-)-0(V)lk~ ^Cst. 

Proof. Set = <p e — <f>. It satisfies 

(8 t + V^ik • V)<4 - \ (|V0| 2 - |V^| 2 ) + \a\ 2 - K| 2 ; w^ t=0 = 0. 

By Proposition 3.1, the right hand side is O(e) in L°° . Integration along the 
characteristics associated to dt + V^ e ik ■ V (see Remark 1.2) yields the result. □ 



The first estimate (1.7) of Theorem 1.3 follows easily: 



r ,„ ' = n (- - a, ' = (a £ - a)e^ /£ + a (e l<t >" ' e - e^/ £ ) 



i<t>/e = a e^4> E /e _ a£ i4>/e = ^ /e 

= 0(e) + ae^'+^/^2ism 



2e 

where the 0(-)'s stand for estimates in L°°([Q,T*\ x K"). 



0(e) + 0(t), 



12 



T. ALAZARD AND R. CARLES 



To improve (1.7) to (1.8), we need the next term in the asymptotic expansion of 
((f> e ,a e ) in terms of powers of e. Introduce the system: 



(3.2) 



d t (f> {1) + V(^eik + <f>) ■ V0 (1) + 2 Re (aa«) = ; </{ t =o = °- 
dta^ + V(0 cik + <j>) • Va« + V0« • Va + A(0 cik + tj>) 



+ \a&4> {1) = '-Aa ; a« = ai . 

It is easy to see that this linear system has a unique classical solution such that 
aW,V^W e C*([0,T»];X s ) for all s > n/2. Reasoning as in the proof of Corol- 
lary 3.3, we see that we have also (j)^ € C([0, T*]; X s ). Moreover, mimicking the 
proofs of Proposition 3.1 and Corollary 3.3, we have the following result, whose 
proof is left out: 

Proposition 3.4. Let (<f) e ,a e ), (</>, a) and (c/)^,a^) be given by (2.6), (3.1) and 
(3.2) respectively. Denote r„ = a E Q — a — ea\. For all s > n/2 + 2, 

||V((/) e -0-e(?!) (1) )|| L oo ([OjTii];Xs) + ||a £ -a-£a (1) || L oo ([O!Tt];Xs) C s (e 2 + \\r E \\x>) ■ 

In addition, there exists C such that if s > n/2 + 2, 

\\(/) £ - 4>- £</> (1) ||L~([o/r»]xR") < C (e 2 + ||rg|| x «) • 

We can now complete the proof of Theorem 1.3: 



= (a £ - a) e^/ £ + a (e^/ £ - ^(W 1 ')/^ 
= 0(e) + ae^+^+^/^sin^ 
= 0(e) + 



2s 



K\\x< 
e 

This yields (1.8), along with Remark 1.5. 

Remark 3.5. Following the same lines, we see that if a§ is known up to order 
0(e N+1 ) in X s for some s > n/2 + 2, JVeN, then we can construct an approximate 
solution v £ N such that 



\\u s -v%\\ L oo {[0iT , ] , xa) =O{e N ). 

To conclude this paragraph, we note that if we know that the initial corrector 
ai is not only in X°° , but in H°° , then Theorem 1.3 becomes more precise. 

Corollary 3.6. Under the same assumptions as in Theorem 1.3, suppose moreover 
that a\ G H°° , and 

|| — a — eai\\ Hs = O (S £ ) , Vs > 0, with S £ = o(e). 

Then (1.8) can &e improved to: 



(3.3) sup u e (t,-)-a(t,-)e** w l*'-)e 

t€[0,T»" 



l 0( 1 '(t,.) pl (0(t,-)+0c, k (t,-))/e 



<5 e 

O 

l°°c\l 2 V e 



Essentially, one just has to notice that the error estimates in Propositions 3.1 
and 3.4 can then be measured in H s instead of X s . Note also that in (3.3), it may 
happen that none of the two functions is in L 2 . 

4. Examples when cik = 
In this paragraph, we consider (1.1)— (1.2), and we assume ^> e ik = 0. 
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4.1. An example from [9]. As an application, we can recover and improve the 
result of [9], in the case of the whole space (the space variable x lies in a bounded 
domain in [9]). Assume that 

og(x) =ao(x) = e iea{x \ $o G H 00 (M. n ; R) ; <jy = V = 0. 

That is, we consider: 

isd t u e + jAu £ = (\u £ \ 2 -l)u £ ; u £ (0,x) = e ie ° {x) . 

Then Oq = a G and we see that: 

• 0=0 and a is independent of time: a(t,x) = a (x) = e ie °( x K 

• <jP-} solves 

dU (1) =Im(oAo), 
so that 6(t,x) := (j>^\t,x) + 8 (x) solves: 

(<9 t 2 -A)6> = ; 6(0,x)=e (x) ; d t 6(0,x)=0. 

Note that (0, a) remains smooth for all time, so we can take T* arbitrarily large 
(see Remark 3.2). Since from Theorem 1.3 and the above corollary, 

sup ||« e (t,-)-a(t,-)e < * <1> ^|Uoc nia = sup \\u e (t, •) - e i6 ^\\ L ~ nL i = 0(e), 
te[o,T„] te[o,T,] 

where > is arbitrary. We recover [9, Theorem 2] in the case of the whole 
space, with no restriction on the space dimension, and a precised error estimate. 
Note also that in view of Remark 3.5, we can justify [9, Proposition 5] (giving the 
e-order corrector for u e ), and get a complete asymptotic expansion for u e . 

4.2. When |ag| 2 — 1 G L 2 . As in Corollary 3.6, assume that (1.6) is precised to 

||oq — a - eai\\ HS = o(e), Vs > 0. 
where do G X°° and a\ G H°°. Assume moreover that 

|a | 2 - 1 G L 2 (R n ). 

Then (2.2) yields: 

j t \\\a £ (t)\ 2 - l\\ 2 L2 = 4 J \\a £ (t,x)\ 2 - l\Re(a £ (t,x)d t a £ (t,x))dx 

< \\\a £ (t)\ 2 - 1\\ L2 \\a £ \\ L ~ (||V$ £ • Wa £ \\ L 2 + ||a £ A$ £ || L 2 + ||Aa £ || i2 ) 

< \\\a £ (t)\ 2 - 1\\ L2 \\a £ f xs (||V$ £ || L ~ + ||A$ £ || i2 + 1) , 
where we consider s > n/2 + 2. Therefore, Proposition 2.1 shows that 

\u £ \ 2 - 1 G C([0,T,];L 2 (K")). 
Note that this property holds even if V — Min ^ 0. 

4.3. When a%(x) ~ 1 as \x\ — > oo. In a spirit similar to [20] (where the authors 
choose 6*o = 0), assume that 1^ = 0, 4>o(x) — v°° ■ x for some v°° G K™, and 



f _ e ie {x) _ £ai 



= 0(S £ ), Vs > 0, where 6 G H°° is real- valued. 



Then as in §4.1, we compute: 

L,oo|2 

0(t, x) = v°° ■ x - ] -^t ; a(t,x) =a a (x-v°°t) = e Uo(x-v°°t) ^ 

We also note that > can be taken arbitrarily large. In addition, we check that 
</>W is such that (t, y) = </>W (t, x + v°°t) solves: 

(d 2 t - A)4> (1) =Im(a Aa ) = A6» ; 4> {1) (0, z) = ; d t 4> {1) (0, a:) = -2 Re (obd) • 
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Therefore, Corollary 3.6 yields 



sup 

te[o,T] 



u £ (t, •) - e 



where 8 is given by 6(t, x) — 0(t, y) 
(d 2 -A)6 = ; 



L°°r\L 2 
with: 



= 0[e + 



y=x—v°°V 

|t=o = #0 ; d t 9\ t=Q — —2 Re (a ai) . 



5. Time propagation of the condition at infinity: (f> cik ^ 

In this section, we assume that |<2,q | 2 — 1 G L 2 (R"), and aim at understanding how 
this condition is propagated on the time interval [0, T»] when O jk =/= 0. Essentially, 
we have \u £ (t, -)| 2 - 1 G i 2 (M") for t G [0, T»] if and only if cik = 0. The function 
</>eik is identically zero if and only if V^ ua d = ^quad = 0: that case was developed in 
§4.2. We compute 

j t ||K(t)| 2 - lf L2 < 4 ||KW| 2 - 1|| L2 ||a £ (t)|| L ~||c^(t)|| L2 . 

In the above estimate, we assumed that dta £ (t, •) G L 2 . Let us now examine this 
condition. In view of Proposition 2.2, we know that all the terms in the second 
equation of (2.6) are in L 2 (M. n ), except possibly d t a £ , V0 c ik • Va e and a £ A(j) c i k . 
Therefore if </> cik = 0, we infer that |u £ (i, -)| 2 - 1 G L 2 {R n ) for all t G [0,T*]. 

Assume now that is not zero. To gather the terms dta £ and V^ e ik ■ Va £ 
together, consider the change of variable of Remark 1.2, and set 

a £ {t,y) = a £ (t,x(t,yj). 

Since the Jacobi determinant det V y x(t, y) > is bounded from above, and from 
below away from zero for t G [0, T*] C [0, T], d t a £ (t, •) and d t a £ (t, •) are simultane- 
ously in i 2 (M"). Given A0 c ;k is a function of time only, we have 

d t a £ = -is e A0 ci k + C([O,T»] ;J L 2 ). 

We are in a case where a e A</> oi k ^ L 2 . To overcome this issue, consider 



2 


2 




- 1 








L 2 





a e( t ) e / 'lVQ(r)dr 



L 2 



where Q is given by Lemma 1.1. For t G [0, T»], this quantity is equivalent to: 



L 2 



We have: 

d_ 

dt 



2 


2 


- 1 






L 2 



a £ (t)eti TrQ( ^ dT 



\a £ (t)\\ L ^x 



L 2 



d t a £ (t) + -a £ (t)A0 cik (t) 



L 2 



\\c infer that o^e-^QW*- - 1 G C([0, T,]; i 2 ), hence 

u e^ e / 'TrQ(r)dr 2 _ j g (7([0, T,]; L 2 ). 

Morally, for i G [0, T*], the modulus of w £ goes to cxp(— J Q Tr Q(r)dr) as |x| — > oo. 
We conclude by some examples that illustrate this analysis. 
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Example 1. Consider the case where qU ad = 0, and V^uad(^) = w 2" ^ s an 
isotropic harmonic potential (ui > 0). Then we compute 



Id 2 

<^ik(t, x) = tan(wt), i e [0, T] c 



2/2 



and cxp ^— y Tr Q(r)c?r^ = cxp ^-y- ^ tan(wr)dr^ = (cos(wt)) n l 2 _ 

Therefore, the "limit of the modulus of u e at infinity" grows at time evolves. If in 
Proposition 2.2, we can take T* arbitrarily close to 7r/(2w), this suggests that there 
is some sort of "blow-up at infinity" at t approaches tt/(2u>). 

Example 2. Consider the case where </> qua d = 0, and V qua d{ x ) = ^ uj2 ^2~ 1S an 
isotropic repulsive harmonic potential (u > 0). We have 

M 2 

</> e ik(i, x) = u—^- tanh(wi), i € [0, +00 [, 
and exp (- Tr Q(r)dr^ = (cosh(wt))"" /2 . 

Therefore, the "limit of the modulus of u e at infinity" decays at time evolves. 
Example 3. Consider the case </> qua d = — M 2 /2, and Kiuad^) = 0. We compute 

<j>eik(t,x) = g^l^ . *e[0,l[, and exp^-^ Tr Q(t)c?t^ =(l-t)~ 
This case is similar to the first example. 

Example 4- Consider the case </> qU ad = M 2 /2, and V qua <i(x) = 0. We have 
<f> cik (t,x) = 2 ^ +1 y te[0,+oo[, and cxp ^- jf Tr Q(r)dr^ = (l+tp l/2 . 
This case is similar to the second example, provided that we consider positive times. 

6. On the hydrodynamic limit 

In this paragraph, we consider the setting of cither Theorem 1.3 or Theorem 1.10. 
That is, the semi-classical limit is justified for small time in Zhidkov spaces. Let 
$ = 4>dk + 0, v = V<I> and p — \a\ 2 . As is easily checked, (p, v) solves the following 
compressible Euler equation: 

, . ( dtp + div(pv) = ; p t =o = \a \ 2 . 

[dtv + v Vv + W + V/(p) =0 ; V| t=o = V0 o , 

where /(/>) = p — 1 in the cubic case, and /(p) = p 2 + Xp in the cubic-quintic case. 
To simplify the discussion, assume in this paragraph that V^ ua( j — 0quad 

= 0, hence 

4>cik = 0. Proposition 3.1 implies in particular the convergence of the main two 
quadratic quantities, as e — > 0: 

• Density: |w e | 2 -> p in L°°([0,T*] x K"). 

• Momentum: Im(ew £ Vu e ) -> pv in L°°([0,T*] x R"). 

It should be noted that if we assume only that for some s > n/2 + 2, 

|| fl o — a o||x» = <^o = as e ^ 0, 

the proof of Proposition 3.1 shows that we have: 

||V(> £ - </>)\\l°°([0,T.];X*) + IK - a||l,°°([0,T„];X«) = O (e + *o) • 

Therefore, 

\u e \ 2 = p + O {e + 5 e Q ) ■ lm(eH e \7u e )=pv + 0(e + S E a ). 
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To have a more precise asymptotics, it is necessary to work with the assumption of 
Theorem 1.3. If for some s > n/2 + 2, 

||oq — a — sai\\x s = 5f = o(e) as e — > 0, 

we get: 

\u e \ 2 =p + 2e Re (aa^ + O (e 2 + 5{) . 
lm{eu e Vu e ) = pv + e (2 Re (aa (1) ) v + pV^ (1) ) + O (s 2 + <Jf) . 

Finally, note that in general, even if a\ = 0, the modulation <j>^ is not trivial. 
Suppose that 01 = 0: (3.2) shows that d t a^L Q 7^ 0, because of the source term 

|Aa. Therefore, even if <f>^L Q = ^t^|t=o = 0' we nave ^^[t=o 7^ in general, and 
the correctors of order s in the above asymptotics are not trivial. 

However, if a is real-valued and a\ = 0, then a is real-valued, a^ 1 ' is purely 
imaginary, so cjP^ = 0. The same holds if a is real-valued and a\ is purely imagi- 
nary. 

We end this section by studying the hydrodynamic limit in the case when SlcK" 
is a regular domain with bounded boundary dfl and n G {2, 3} (either a bounded 
domain or an exterior domain). To simplify the presentation, we consider the case 
without external potential and without linear or quadratic initial phase. The Gross- 
Pitaevskii equation is then supplemented with the Neumann boundary condition: 



(6.2) 



F 2 

ied t u E + —Au £ = (\u e \ 2 -l)u £ in ft, 

— — =0 on 9ft, 

on 



where n is the unit outward normal to 9ft. Consider the corresponding limit system 

!d t p + div(pV^) =0 in ft, 

3 t 0+^|V0| 2 +p-l = O in ft, 
V</> • n = on dn. 

In [20], Lin and Zhang proved that if n = 2, then the quadratic observables |u £ | 2 
and elm(u e Vu e ) converge towards the density p and the momentum pVcj). In the 
spirit of the pionnering work of Brenier [6], in [20] the strategy of the proof is to 
estimate the modulated energy functional 



E e :=\ f \eVu e - m £ V0| 2 + (\u e \ 2 - p) 



2 dx. 



The assumption n — 2 does not enter into the analysis of E e and only corresponds 
to the fact that they used the Brezis-Gallouet inequality (see also [9]) to define 
sufficiently smooth solutions to the Gross-Pitaevskii equation. They are now several 
3D results (see [16, 4, 11, 12]), and hence one can justify the hydronamic limit for 
n e {2,3}. In particular, Theorem 6.1 below is not new, but rather an update. 
Yet, our main purpose here is to establish a local version of the modulated energy 
functional. This is done in the proof of Theorem 6.1 (see (6.4)), by following the 
approach introduced in [3]. 
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Theorem 6.1. Letu £ and (p,(f>) be classical solutions of (6.2) and (6.3) satisfying, 
for some fixed T > 0, 

u £ e C([0,T]:X 2 (Q)), \u £ \ 2 - 1 e C([0,T];L 2 (n)), 

P eC([0,T]:X\Q)), p-le C([0,T];L 2 (ft)), 

V<t>, V 2 </>, V 3 e C([0, T]; L 2 (f2) n L°°(fi)). 

Assume that initially 

\\eVu e - iu§V0o|| ia(n) + ||K| 2 - Po|| i2(a) - 0(e), 

i/ien 

K| 2 -p = 0( £ ) in L°°([0,T];L 2 (r!)), 
eIm(u e Vu £ ) -pV0 = O(e) in L°°([0, T]; Lj, c (fi)). 
Proof. The idea consists in filtering out the oscillations by the change of unknown 

a £ (t,:r) := u £ (t, aOe"^*'^. 

The amplitude a £ solves 

<9 t a £ + • Va £ + \a £ A(j) - ^Aa 6 = -- (|a £ | 2 - p) a 6 . 

Next set 
We easily find that 

d t q £ + div (Im (a £ Va £ )) + div(<fV</>) = 0. 
Furthermore, with this notation, the equations for ip £ := Va E read 

d t ip £ + V<£ • VV> £ + ^ £ Acj) + ip £ • VV0 + ia £ VA0 + ^ £ i/> £ + ia £ \7q £ = i S -^ £ . 

Also, note that ip £ ■ n = e"^ £ (V« E - ie- 1 u £ V(f>) ■ n = on dtt. 
We now introduce the modulated energy 

e £ := |^ £ | 2 + (<f) 2 . 

The key point is that 

g £ div(lm(a £ Va £ )) + Re (ia £ (Vq £ ) ■ = div ( lm(q £ a £ xp £ )) . 
Hence, directly from the previous equations, we have 

(6.4) d t e £ + div(e £ V0) + div(2 lm(q £ a £ ^ £ )) + div (e Im (ijf • V?/> £ ) ) 

= -(q e ) 2 A<P - Re ((2t/> £ • VV</» + a £ VA</>) • ^) . 

We claim that 

(6.5) E £ (t) = ||e £ (t)|| il(0) < ||e £ (0)|Ui ( o) cxp (Ct) + C, 

for some constant C independent of e. Since v ■ n = and tp s ■ n — on dQ, 
by integrating in space and using the Gronwall's lemma, to prove (6.6), the only 
delicate point is to prove that, 

(6.6) J a £ VA0-^ e dx ^C\\e e \\ L i m +C. 

To do so, as in Lemma 1 in [16], let x € Co(C) be such that < \ < 1, x( z ) = 1 
for \z\ < 2, and \{z) — for \z\ ^ 3. Then write a £ — b £ + c £ where 6 £ = x(a £ )a £ 
and c £ = (1 - x{a e ))a e - We have |6 £ | 3, |c £ | < ||a £ | 2 - l| and hence 

W\\ L -<fl) O, ||c £ || L2(0) < |||a £ | 2 -l|| L2(f2) <e|k e || L2( o) + llp-l|lL 2 w 
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The desired estimate (6.6) then follows from 



b E VA(j) ■ </> 
c £ VA0 • ~$ 



<||6 £ || iMW ||VA0|| i2(fi) ||^|| i2 



L!(f2) 



< c £ 



L 2 (n) 



||VA0|| 



L°°(fi) 



L 2 (Q) 



and the elementary inequality -^/a; $C 1 + x. 
Since 



e £ = — \eVu £ ~iu £ V<j)\ 2 + — (\u £ \ 2 - p) , 

the family (e e (0)) eg ] ^ is bounded in by assumption. Consequently, it fol- 

lows from (6.5) that (e £ ) Ee]al] is bounded in L°°([0, T]; L 1 (0)). 

By definition, this implies that |u e | 2 -p = 0(e) in L°°([0, T];L 2 (f2)). It remains 
to prove that 



Write 



elm(u £ Vu £ ) - pV(j> = 0(e) in L°°([0, T\; Ly oc {Sl)) . 



eIm(u £ Vii £ ) - pV0 = eIm(a £ Va £ ) + (|a £ | 2 - p) 



Since V</> € L°°([0,T] x f2), the previous result implies that the second term is 0(e) 
in L°°([0, T]; L 2 (£l)). With regards to the first one, again write a £ = b £ + c £ and 
use the obvious estimates 



||£lm(6 e Va e )|| I ,2( n) ^ e||6 £ ||i,c»(a)||Va £ || z ,2 (n) ==S 3e ||e £ 



1 1/2 
1/2 



||eIm(c £ Va £ )|| il( o) < e\\c £ \\ LHn) \\Va £ \\ L i (n) ^ Ce\\e £ \\^ {n) 
This completes the proof. 



-Ce 2 \\e £ 



□ 



7. CUBIC-QUINTIC NONLINEARITY 

In view of Theorem 1.10, we now consider (1.9) in the case where the elliptic 
region becomes relevant: A < 0, and assume for instance that there exists x £ M. n 
such that |ao(x)| 2 < |A|/2. If we write u £ — a £ e l * e/,£ , where (a £ ,$ £ ) is given by 
(2.2), then we naturally have to consider the limit system: 



(7.1) 



\m\ 2 



fx 







d t a + Mcj) ■ Ma + -aA(j> = 



\t=o 



it=o 



Setting v = V<fi, we find: 

dtv + v Vv + Vfx (\a\ 2 ) =0 



(7.2) 



d t a + v ■ Va + -adivv = 



\t=o 



a . 



a . 



Then [22, Theorem 3.2] shows that (7.2) is strongly ill-posed in Sobolev spaces. 
The problem remains in Zhidkov spaces, since analyticity is essentially necessary. 
Indeed, Hadamard's argument (sec [22] and references therein) shows for instance 
that if (j>o is analytic near x, then (7.2) has a C 1 -solution only if ao is also analytic 
near x. So it may happen that (7.2) has no solution in X s , even for s large. 

On the other hand, if one is ready to work with analytic regularity, then it 
becomes possible to justify the semi-classical limit for (1.9); see [15, 27]. 
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